B = 3 Tetrahedrally Symmetric Solitons in the Chiral Quark Soliton Model 
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In this paper, B — 3 soliton solutions with tetrahedral symmetry are obtained numerically in the 
chiral quark soliton model using the rational map ansatz. The solution exhibits a triply degenerate 
bound spectrum of the quark orbits in the background of tetrahedrally symmetric pion field con- 
figuration. The corresponding baryon density is tetrahedral in shape. Our numerical technique is 
independent on the baryon number and its application to B > 4 is straightforward. 
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The Chiral Quark Soliton Model (CQSM) was devel- 
oped in 1980's as an effective theory of QCD interpolat- 
ing between the Constituent Quark Model and Skyrme 
Model [0, 0|. In the large N c limit, these models are 
identical 0^ 

The CQSM is derived from the instanton liquid 
model of the QCD vacuum and incorporates the non- 
perturbative feature of the low-energy QCD, spontaneous 
chiral symmetry breaking. The vacuum functional is de- 
fined by; 



Z = \ VirVipVip 1 exp 



i I d xijj (i</) — mlf 75 ) ip 



(1) 



where the SU(2) matrix 



ins = -^y^-u + ^—^-u^ with u = cx p (* T? ' 7? //^) 

describes chiral fields, ip is quark fields and m is the dy- 
namical quark mass. f n is the pion decay constant and 
experimentally f n ~ 93McV. Since our concern is the 
tree-level pions and one-loop quarks according to the 
Hartree mean field approach, the kinetic term of the 



J 



pion fields which gives a contribution to higher loops 
can be neglected. Due to the interaction between the 
valence quarks and the Dirac sea, soliton solutions ap- 
pear as bound states of quarks in the background of self- 
consistent mean chiral field. N c valence quarks fill the 
each bound state to form a baryon. The baryon number 
is thus identified with the number of bound states filled 
by the valence quarks [f| . 

For B = 1 and 2, the spherically symmetric soliton 
H 0, and the axially symmetric soliton || were found 
respectively. Upon quantization, the intermediate states 
of nucleon and deuteron between the Constituent Quark 
Model and Skyrme Model were obtained. 

The vacuum functional in Eq. (Jl|) can be integrated over 
the quark fields to obtain the effective action 



S e e[U] = -iN c \ndet(i$-mlI r 
i 



(2) 
(3) 



= --iV c SplnD t L» 

where D = ifi — mil 15 . This determinant is ultraviolet 
divergent and must be regularized. Using the proper- 
time regularization scheme, we can write 



= >c 



'l/A 2 



-D*Dt 



= 2^ T 



2^ 



'l/A 2 T 



e -r(H^ 



(4) 



where T is the Euclidean time separation, A is a cut-off 
parameter evaluated by the condition that the derivative 
expansion of Eq. (|^) reproduces the pion kinetic term with 
the correct coefficient i.e. 

/J = f#f^-\ (5) 
4tH J 1/A 2 t 

and H is the Dirac one-quark Hamiltonian defined by 
a-V 

H = —— + (3mU^ . (6) 



D = D(U = 1) and H = H(U = 1) correspond to the 
vacuum sectors. 

At T — * oo, we have e lSoff ~ e _l£;scaT . Integrating over 
ui in (0) and constructing a complete set of eigenstates of 
H with 

H]v) = e v ]v) , H Q \v)W=<M\v)W , (7) 
one obtains the sea quark energy B 
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1 f°° 
E sea [U] = -r-=N c / 



Z3J2 



(8) 



In the Hartree picture, the baryon states are the quarks occupying all negative Dirac sea and valence levels. Hence, 
if we define the total soliton energy E to tai, the valence quark energy should be added; 



E total [U] =N C J2 [U] + E sea [U] 



(9) 



where E l val is the valence quark contribution to the i th baryon. 

The baryon density (60) for the baryon number B soliton is defined by the zeroth component of the baryon current 



(bo) 



N r B 



N r B 



£(n^(e„)+ S ign(e„)Af(e„))<H^^ 

V V 

I 



(10) 



where 



1 



-i/U^A 



A/"(e„) = 



and n v is the valence quark occupation number. 

For B = 3, it is expected the solution to have a tetra- 
hedral symmetry from the study of the Skyrme model 
|UJ . Therefore, we shall impose the same symmetry on 
the chiral fields using the rational map ansatz. According 
to the ansatz, the chiral field can be expressed as Q| 



where 



nR 



l + \R(z) 



U(r,z) = exp(if(r)n R -r) 



(2Re[i?(z)],2Im[i?(z)],l- \R{z)\ l 



(11) 



and R(z) is a rational map. 

Rational maps are maps from CP{\) to CP{\) (equiv- 
alently, from S 2 to S 2 ) classified by winding number. It 
was shown in |ll[] that B = N skyrmions can be well- 
approximated by rational maps with winding number 
N. The rational map with winding number N possesses 
(2./V + 1) complex parameters whose values can be deter- 
mined by imposing the symmetry of the skyrmion. Thus, 



B — 3 rational map with tetrahedral symmetry takes the 
form; 



R(z) = 



z{z 2 - VSi) 



(12) 



where the complex coordinate z on CP(1) is identified 
with the polar coordinates on S 2 by z — tan(0/2)e I¥ ' 
via stereographic projection. Substituting ( |l2| ) into (pi]), 
one obtains the complete form of the chiral fields with 
tetrahedral symmetry and winding number 3. 

Apparently, the chiral fields in ([ll]) takes a spherically 
symmetric form. Therefore one can apply the numeri- 
cal technique developed for B = 1 to find B = 3 with 
tetrahedral symmetry 0. 

Demanding that the total energy in ([)]) be stationary 
with respect to variation of the profile function /(r), 



6f(r 



■E, 



total 



, 



yeilds the field equation 



S(r) sin f(r) — P(r) cos/(r) 



(13) 



where 



S(r)=N c J2{^(^) + ^(euW(eu))(v\l°S(\x\-r)\ V ) (14) 

P(r) =7V c ^(n^(e,)+sign(e,)AA(e,))(^|i 7 7 5 ^-^(kl -r)\v) . (15) 

I 

The procedure to obtain the self-consistent solution of Eq.(|l3"|) is that 1) solve the eigenequation in (p]) under 
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FIG. 1: Spectrum of the quark orbits are illustrated as a function of the "soliton size" X. Trial profile function is defined 
by f(r) = -7T + Tvr/X for X < r and f(r) = for X > r. The state coupled with (K, M) = (0,0) is denoted as S and 
(K,M) = (1, M) as Pm {M — ±1,0) with parity ±. The energy of the three P states are all degenerate. 



TABLE I: A schematic picture of the matrix elements 
A(K'M',KM) up to K, K' = 2. 5", Pi, P and P_i refer 
to the elements coupled with (K, M) = (0,0), (1,1), (1,0) 
and (1, —1) respectively. Other elements are all 0. 



(0 0) (1 1) (1 0) (1 -1) (2 2) (2 1) (2 0) (2 -1) (2 -2) 



(0 0) s s 

(1 1) Pi Pi 

(10) P P P 

(1 -1) P-i P-i 

(2 2) Po P Po 

(2 1) P_i P_i 

(2 0) 5 S 

(2 -1) Pi Pi 

(2-2) Po Po Po 



an assumed initial profile function fo(r), 2) use the resul- 
tant eigenfunctions and eigenvalues to calculate S(r) and 
P(r), 3) solve Eq.([l3|) to obtain a new profile function, 
4) repeat 1) — 3) until the self-consistency is attained. 

For convenience, we shall take fa(r) = — ne~ r / x . To 
solve Eq.(0), we construct the trial function using the 
Kahana-Ripka basis Q ; 

mam i=l K=0 M=-K 

where \E' + and ^~ stand for parity (— 1) and (— 1) K+1 



respectively, tp is the Kahana-Ripka basis and K is the 
grand spin operator which is a good quantum number 
in the case of B = 1 hedgehog. The basis is discretized 
by imposing an appropriate boundary condition for the 
radial wavefunctions at the radius r max chosen to be suf- 
ficiently larger than the soliton size. And also, the basis 
is made finite by including only those states with the mo- 
mentum k as k < k m ax- The results should be, however, 
independent on r max and k max . 

According to the Rayleigh-Ritz variational method 
|l2| , the upper bound of the spectrum can be obtained 
from the secular equation for each parity; 

det (A± - eB ± ) = (17) 

where 

A ± (K'M / , KM) = J2 I dW^M-H^M 
ij=i J 

K±{K>M>,KM) = E / d3 ^K'W^M ■ 

y=l J 

For K — > co, the spectrum e becomes exact. Eq.(|l7|) can 
be solved numerically. 

Since the chiral field in Eq.(pT|) is less symmetric than 
the B = 1 hedgehog, the hamiltonian has no grand spin 
symmetry. As a result, the states with different grand 
spin couple strongly and level splittings within the K 
blocks occur. In Table 1 we present the schematic picture 
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FIG. 2: Surface of the baryon-number density with bo = 
0.4 fm -3 . 




r[fm] 



FIG. 3: Profile function f(r) of the rational map ansatz for 
5 = 3, and of the hedgehog ansatz for B = 1. 



of the matrix elements A(K'M', KM). Although the size 
of the matrix A(K'M', KM) becomes quite large, due to 
the symmetry of the chiral fields, the functional space can 
be rearranged to reduce the size. Consequently, the space 
is divided with four blocks for each parity. 

Fig. |l| shows the spectrum of the quark orbits as a 
function of the soliton-size parameter X . The P + orbit 
diving into the negative energy region is triply degener- 



ate. As discussed in baryon number of the soliton 
equals to the number of diving levels occupied by N c va- 
lence quarks. Thus putting N c — 3 valence quarks on 
each of the degenerate levels, one obtains the B = 3 soli- 
ton solution. 

Fig. ^| shows the corresponding baryon density. As can 
be seen, it is tetrahedral in shape. Therefore, we confirm 
that the lowest lying B — 3 configuration is tetrahedrally 
symmetric. This result is consistent with the B = 3 
skyrmion obtained by Braaten et. al JlO| . 

Fig. |^ shows the self-consistent profile function. For 
the total energy of the solution we obtain E to tai = 3596 
MeV which is almost comparable to three times of the 
B = 1 mass, and root mean square radius is \J (r 2 ) ~ 0.6 
fm. Our soliton seems to be tight object. This is mainly 
due to the missing of higher components of K in our 
calculation. Their contribution becomes significant near 
the surface of the soliton and hence inclusion of the higher 
components will improve the size of the soliton. 

Finally, we would like to mention that our result ver- 
ifies the validity of the rational map ansatz for the Chi- 
ral Quark Soliton Model. The numerical technique used 
here is quite general and its application to B > 4 will be 
straightforward. 
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